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13  ABSTRACT 

This  paper  is  concerned  with  the  investigation  of  the  inelastic 
buckling  of  a  deep  spherical  shell  subject  to  a  uniformly  distributed 
external  pressure.  The  geometry  of  the  shell  is  considered  to  be  axisymmet-j 
rical  while  the  shell  thickness  may  vary  as  a  function  of  the  polar  angle. 
The  edge  of  the  shell  is  supported  elastically.  The  material  of  the  shell 
is  assumed  to  satisfy  the  generalized  Ramberg-Osgood  stress-strain  relations 
and  a  power  law  of  steady  creep.  The  analysis  is  based  on  Sanders’  non¬ 
linear  theory  of  thin  shells  expressed  in  an  incremental  form  and  Hill's 
theory  of  inelastic  bifurcation.  Computations  are  carried  out  by  a 
numerical  iterative  procedure  associated  with  a  finite  difference  method. 
Solutions  are  sought  for  both  the  axisymmetrical  inelastic  buckling  and 
the  asymmetrical  bifurcation. 
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ABSTRACT 

This  paper  is  concerned  with  the  investigation  of  the  inelastic 
buckling  of  a  deep  spherical  shell  subject  to  a  uniformly  distribu¬ 
ted  external  pressure.  The  geometry  of  the  shell  is  considered  to 
be  axisymmetrical  while  the  shell  thickness  may  vary  as  a  function 
of  the  polar  angle.  The  edge  of  the  shell  is  supported  elastically. 
The  material  of  the  shell  is  assumed  to  satisfy  the  generalized 
Ramberg-Osgood  stress-strain  relations  and  a  power  law  of  steady 
creep.  The  analysis  is  based  on  Sanders’  nonlinear  theory  of  thin 
shells  expressed  in  an  incremental  form  and  Hill’s  theory  of  ine¬ 
lastic  bifurcation.  Computations  are  carried  out  by  a  numerical 
iterative  procedure  associated  with  a  finite  difference  method. 
Solutions  are  sought  for  both  the  axisymmetrical  inelastic  buckling 
and  the  asymmetrical  bifurcation. 
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1.  INTRODUCTION 

In  recent  years.,  new  developments  in  deep  ocean  exploration 
demand  further  refinement  and  sophistication  in  the  design  of  sub- 
mersibles.  An  important  problem  encountered  in  the  design  of 
submerged  hulls  is  the  selection  of  a  suitable  material  with  high- 
strength  and  low-density  properties.  In  view  of  the  nature  of  the 
material  available  and  due  to  the  limitations  on  structural  weight, 
it  is  necessary  that  an  accurate  stress  analysis  of  the  hull  and 
a  reliable  prediction  of  the  collapse  load  be  established. 

When  the  submerged  depth  is  large,  the  configuration  of  the 
hull  is  usually  spherical,  fhe  spnerical  shell  can  provide  a  good 
structural  layout  from  the  viewpoint  of  stress  analysis.  Further¬ 
more,  it  has  the  advantage  of  low  drag  during  the  motion  of  the 
hull.  CNR’s  vehicle  ALVIN  is  a  typical  example  of  the  spherical 
hull.  For  this  type  of  pressure  hull,  when  it  is  submerged  to 
a  large  depth,  the  deformation  of  the  hull  may  become  finite  and 
the  state  of  stress  in  the  hull  may  reach  the  inelastic  range. 

In  a  critical  condition,  structural  failure  may  occur  as  a  result 
of  inelastic  buckling. 

Since  Shanley  first  introduced  the  concept  of  inelastic  bi¬ 
furcation  f  c  t  column  under  increasing  axial  compressive  load  [1], 
substantial  T  rjgress  has  been  made  in  the  field  of  inelastic  buck¬ 
ling  of  structures.  Research  along  this  line  can  be  found  in  a  com¬ 
prehensive  survey  paper  prepared  by  Sewell  [2].  The  general  study 
of  the  inelastic  stability  of  solids  based  on  the  concept  of 


uniqueness  of  solution  was  given  by  Hill  [3,4].  He  concludes  that 
Shanley’s  concept  of  inelastic  buckling  under  increasing  load  is 
a  natural  sequel  to  his  general  theory.  In  Hill’s  work,  a  varia¬ 
tional  principle  is  established  for  finite  deformation  of  inelastic 
solids  based  on  the  existence  of  convex  plastic  potential  surface. 
This  variational  principle  is  essentially  an  extension  of  the 
well-known  principle  of  minimum  potential  energy  in  finite  elas¬ 
ticity.  Hill’s  variational  principle  incorporated  with  Rayleigh- 
Ritz  technique  or  Galerkin  technique  has  been  employed  for  the 

analysis  of  inealstic  buckling  of  thin  shells  [5,6,7]. 

» 

Similar  tc  the  case  of  elastic  buckling,  the  critical  con¬ 
dition  for  inelastic  buckling  is  usually  sensitive  to  the  initial 
imperfections  existing  in  the  structural  geometry  or  loading  con¬ 
dition.  In  fact,  for  many  structures,  such  as  columns  and  complete 
spherical  shells,  the  inelastic  bifurcation  occurs  only  when  the 
geometry  of  the  structure  and  the  loading  condition  are  perfect. 
With  imperfections,  bifurcation  disappears  although  the  collapse 
load  can  still  exist.  It  is  found  that  the  magnitude  of  the 
collapse  load  depends  strongly  on  the  imperfection  of  the  struc¬ 
ture  [8,9],  This  behavior  of  imperfection-sensitivity  indicates 
that  the  analytical  results  derived  from  the  variational  principle 
may  depend  on  the  assumed  deflection  function  used  in  the  Rayleigh- 
Ritz  technique.  To  remedy  the  deficiency  introduced  by  the 
Rayleigh-Fitz  technique,  the  problem  of  inelastic  stability  of 
shells  can  be  investigated  directly  from  the  governing  nonlinear 


shell  equations.  Typical  examples  of  this  approach  are  found  in 

[9,10]. 

In  this  paper,  we  shall  deal  with  the  problems  of  axisymmetri- 
cal  collapse  and  asymmetrical  bifurcation  of  a  deep  spherical  shell 
in  the  inelastic  range  introduced  by  a  uniform  external  pressure. 

s 

The  model  of  spherical  shell  use  in  our  analysis  is  chosen  to  be 

I 

close  to  the  actual  shell  configuration  of  the  ALVIN  vehicle.  The 
thickness  distribution  of  the  shell  is  regarded  as  axisymmetrical 
but  variable  as  a  function  of  the  polar  c.  The  edge  of  the 

shell  is  considered  to  bo  elastically  supported.  The  shell  material 
is  assumed  to  be  isotropic  satisfying  Ramberg-Osgood  stress-strain 
relations  and  a  power  law  of  steady  creep.  Our  analysis  is  based 

T 

on  the  Kirchhoff  assumption  with  small  strains.  Sanders  nonlinear 
theory  for  finite  deformation  of  thin  shells  [11]  will  be  employed 
in  our  formulation. 

2 .  BASIC  EQUATIONS 

Consider  a  deep  spherical  shell  with  radius  of  the  middle 
surface  R  subjected  to  a  uniform  external  pressure  P.  Let  us 
introduce  surfs ce  polar  spherical  coordinates  with  polar  angle  9 
and  longiucde  0  as  shown  in  Figure  1.  In  the  following,  we 
shall  take  6  ar.^  0  at  the  parametric  coordinates.  Thus  our 
lines  of  curvature  are  meridional  lines  and  circles  perpendicular 
to  the  polar  axis.  The  edge  of  the  shell  is  specified  by  the 
circle  0  =  9q.  We  shall  assume  that  the  thickness  distribution 
of  the  shell  is  axisymmetrical,  i.e.  h  --  h(0). 
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Let  us  denote  the  components  of  displacement  in  the  meridional, 
circumferential  and  outward  normal  directions  by  U^,  U2,  and  W 
respectively.  According  to  Sanders’  nonlinear  theory  of  shells 
[11],  the  components  of  rotation  in  the  directions  of  tangents  and 
normal  to  the  middle  surface  are  given  respectively  by 


*i  -  i<ui  -  V’ 


cp2  =  jjr  (U2  -  W^2  esc  0), 

'P  =  2E  (U2  Got  9  +  U2,l  “  Ul,2  CSC  9)' 
where  (  )  ^  =  |~  (  )  and  (  )  2  B  (  )•  The  membrane  strains 


(1) 

(2) 

(3) 


are 


E11  =  K  (W  +  Ul,l}  +  \  (cpl  + 


(4) 


E22  =  £  (W  +  Ui  cot  0  +  U2  2  esc  9)  +  (cp2  +  $"),  (5) 


E12  =  ^R^2  l  +  ^l  2  CSC  9  ”  c°t  9  +  ^  C0tCDo).  (6) 


12' 


The  bending  strains  are 


Kli  “  E  cpl,  1  * 

n 

^22  ~  E  ^2  2  8  cot  9), 


K12  ~  ^®2,1  +  'Pi, 2  osc  0  -  3?2  cot  9)* 


(7) 

(8) 
(9) 


Let  us  denote  the  components  of  membrane  force,  transverse  shear 
and  moment  resultant  by  Nn,  N22,  N12,  Q1,  Q2,  Mn,  M22  and  M12- 
The  equations  of  equilibrium  of  the  deformed  shell  are 


nf-wiPA.  a 


■-  a iMkdtm tlU.  -  ■  ■■  ■  * 


*  ‘ 
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N11  1  +  N11  Cot  9  +  N12  2  CSC  9  "  N22  COt  9  +  'V^°1N11+  co2N12^ 

-  \  Cco(N11  +  N22)]  2  CSC  9"RP  ®l  =  0  *  (10) 


N22,2  +  2N12COS  0^*12  2 sin6  +  ^2Sin  e"^cp2N22  +  cplN12)Sin0 


+  \  [®<N11-H^22)]3  isin0-R?  qv,  sin  e  =  C, 


^1,1^1  COt  6  +  ^2,2  CSc6"(Nll+N22)_(cplNll  +  ^L-2^,1 


“(cplNll'!'cp2N12^cote  “  (pPiNi2  +cp2i'r2?.), 2  esc  6  “  Rp  =  °> 


111  1+  Mncot  6  +  M12,2CS0  9  “M22  COt  6  “  =  °> 


M22  2  +  2M12  COS  9  +  M12  1  Sin  sin  0  =  °* 


These  equations  are  self-consistant  in  a  sense  that  the  principle 


of  virtual  work  can  be  satisfied.  Note  that  nonlinearity  is 


introduced  through  the  strain-displacement  relations  and  equations 


of  equilibrium. 


Ey  the  Kirchhoff  assumption,  the  strain  components  E. .  (i  =  1,2; 

13 


j  =  1,2)  at  any  point  in  the  shell  with  a  distance  z  measured  from 


the  middle  surface  can  be  expressed  in  terms  of  the  membrane  strains 


and  bending  strains  as 


E. .  =  E.  .  +  z  K.  .  , 
13  13  13 


In  Eq.  (15),  we  consider  that  the  shell  is  sufficiently  thin. 


Hence  only  the  linear  term  in  z  is  retained.  Denote  the  stress 


T 


r  .  *-.*>** 
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components  at  any  point  in  the  shell  by  cr^Ci  =  1*2;  j  -  1,2). 

In  the  following  we  shall  assume  that  the  strains  are  small  although 
the  deformation  of  the  shell  may  be  finite.  The  constitutive  re¬ 
lations  can  be  written  as 


“*  • 

- 

- 

» 

E11 

C11  C12  C13 

all 

D,  ! 

X 

E22 

_ 

ore 

U21  u22  u23 

°22 

+ 

D2 

i 

» 

i 

E12 

c  C  c 

C31  ~32  u33 

a12 

o  ! 

..  _ 

—  _j 

- 

J 

where 

C11  =  E  +  I  (2all  “  °22)2’ 

n  1  .  G  ,2 

C22  =  E  +  7  v  2<722  '  °11;  ’ 

p  -  4.  OCr r2 

C33  E“  2Ga12  , 

G 

C12  =  C21  =  ~  E  +  S  ^2all  "  a22^2a22"  all^ 

Q 

C13  =  C31  =  E  (2all"  a22'  °12’ 

G 

C23  =  C32  ®  1  (2ct22  "  all)a12’ 

D1  =  F^ae^all  2  °22^  ’ 

D2  =  F^0e^022  "  2  all^  ’ 


(16) 


(17) 

(18) 

(19) 

(20) 

(21) 

(22) 

(23) 

(24) 


E  is  the  Young's  modulus,  v  is  the  Poisson’s  ratio  and  G  is 
defined  as 

w2(v  E>  for  j2  *  0  a"d  J2  -  <J2W 

0  for  J2  <  0  or  J2  <  (J2)max, 


(25) 
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Et  is  the  tangent  modulus,  is  given  as 


2  2  2 
J2  =  all  +  a22  ~  all  J22  +  3cr12  ’ 


(26) 


(J0L,V  is  the  maximum  value  of  J0  in  its  history,  a  is  the 

^  maX  4  0 

effective  stress  defined  by 


■h 


ae  =  J2  * 


(27) 


The  function  F(ae)  is  given  by 


1  °e  m_1 
F(a  )  =  ±  (-£) 

6  ac  ac 


(28) 


where  m  and  ac  are  material  constants.  The  constitutive 


equation,  Eq.  (16),  is  obtained  by  a  generalization  of  the  uniaxial 
stress-strain  relation: 


m 


v  +  (2  ) 

!=,  {  c  0= 

l  •  m 

i  + 


for  a  a  s  0  ; 
for  a  a  <  0  , 


(29) 


with  the  assumption  that  the  plastic  and  creep  deformations  are 
incompressible.  In  our  study,  we  shall  assume  that  the  material 
of  the  shell  satisfies  Ramberg-Osgood  stress-strain  relation  for 
elastic-plastic  deformation.  Hence  the  tangent  modulus  can  be 
expressed  as 


(30) 


where  oQ  is  the  stress  determined  by  the  intercept  of  the 


uniaxial  stress-strain  curve  with  a  straight  line  through  the 


jgZ2saacsiB^Ms££Mai 


— -  i  i  n  i 
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origin  with  a  slope  0.7E  and  n  is  a  material  constant. 

When  the  shell  is  sufficiently  thin,  the  membrane  forces 


and  moments  are  related  to  the  stress  components  by 

h/2 


Nij  “  5  CTij  dz 


(31) 


and 


-h/2 

h/0 


M.  .  =  C 
13  J 


Z(j^ .  dz  . 


(32) 


-h/2 


Consider  hQ  a*  a  reference  thickness.  Let  us  introduce  the 


following  dimensionless  quantities: 


U, 


'I  W 

U,  —  j  w  = 


R 


'i-TT 

h 
o 


e—  =  TL  Eij'  Kij  "  ^ij  '  (33) 


o  o 


N.  . 


M. 


X=! r  •’T1=Hn-’  nij=^  ’  mii=“T'  (34) 
o  Ehr,  Ehr 


‘o 


„  \ 


13  ±j'  13 


n  —  £»  —  r  ~  12  (•3c:\ 

'  P  ~  EH  '  eii  “  TL  Eii'  Tij  “  TT>  (35) 
0 


°c  .  °o  „  Et 

T"  =  r~  >  T-  =  v-  >  e-  *  v-  > 


ro  -  E~ 


(36) 


Our  governing  equations  for  the  general  deformation  of  the  shell 
can  be  expressed  in  terms  of  these  dimensionless  quantities  as 


®i  =  xc«i 

-  w?1)  , 

(37) 

ec^  =  X(  ^2 

-  W  ?  CSC  0)  , 

(38) 

®  =  ^  ( 1^2 

cot  0  +  U2^x  -  Uj^2  CSC  0)  , 

(39) 

‘'ll  "  "  + 

Ul,l  +  ?X(®1  +  ®2)  ’ 

(40) 

-  j 


UND- 73-7 

9 


-  9  - 


1  2  2 

e22  =  w  +  ui  cot  6  +  u2  2  CSC  0+  ?  0  ^ 


e12  =  -j(u2  1  +  u,  2  esc  e  -  u2  cot  e  +  ~  cp-jcp2). 


K11  ~  *1,1  ’ 


K22  =  ®2,2  CSC  0  +  al  COt  0  J 


(41) 

(42) 

(43) 

(44) 


Kt  rt  — 


1 


v12 


(cp2  1  ^  1  2  9  "*  CDg  0 )  > 


(45) 


nH,l  +  nH  cot  9  +  n12,2  CSC  9  -  n22  COt  0  +  ql 

"(®Inll  +  ®2n12)  "  \  Mnn  +  n22^,2  CSC  9  "  P®1  's  °>  (46) 


n22  2  +  2n12  cos  9  +  n12  L  sin  9  +  <\2  sin  9 


-Ccp2  n22  +  <Pini2^  sin  9  +  ?jr  Ccp(ni;L  +  n22^Jlsin  9"  Pcd2  Sin0  “  °5 

(47) 


ll,l  +  °-l  COt  9  +  ^2,2  CSC  9  "  (nll  +  n22> 


"(cDl  nll  +  *2  n,2^1  ”(c?lnll  +  ©2n12)cot  9"(cPlnl2  +  C02n22^2CSC9 


p  =  0, 


(43) 


mll,l  +  mll  COt  9  +  m12,2  CSC  9  "  m22  COt  9  “  x  ql  ^  °  ^ 


(49) 
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m22,2  +  ^Y>  C0S 

e  +  m12ji  sin  9  - 

i  q2  sin 

*N 

o 

II 

(50) 

e, .  =  e. . 
13  13 

+  C 

3  * 

(51) 

r* . 

~ 

r  . 

r  “ 

i  en 

i 

C11 

C12  C13 

T11 

i'1 

i  . 

i  22 

f 

s= 

C21 

C22  C23 

T22 

+  i 

i 

9 

(52) 

. 

el2 

_  _i 

C31 

C32  C33 

*12 

1  o  . 

L  _j 

C11  = 

1 

X 

Cl+I 

(2'th  -  T 

22 

>2]  , 

(53) 

li 

CM 

CM 

C> 

1 

\ 

[1  +  | 

(2r22  -  tu 

)2]  , 

(54) 

C33  « 

1 

X 

(1  +  V 

+  29t12) 

9 

(55) 

C12  = 

[-V  +  f  (2tu  -  t22)(2t22  - 

Tll)]  * 

(56) 

C13  " 

C31  “  X  f  (2t11  " 

T22)  t12 

9 

(57) 

C23  = 

°32  "  X 

|  (2t22  - 

Tll)  T12 

9 

(58) 

k  F<Te)(2Tll  -  T22>> 


F(t.)  (2t00  -  t,-,). 


7k  iV,ey  vt,22  ~  '11; 


2.2  ,  -  2 

T11  "r  t22  "  T11  t22  +  3t12  * 


,4f-  (|  -1) 

g  =  { ,32  et 
0 


for  j2  4  0  and  j2  »  (j2)max  ; 
for  j2  <  0  or  J2  <  (j^  , 


(59) 

(60) 
(61) 


(62) 


- - -»“■  - - - - 
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.h 

Te  ~  ^2  s 

(63) 

l  t  m-1 

F<V  -  i  (-A 

(64) 

c  c 

r  t_  n  -1-,  -1 

at  =[1  +  7  V?f>  ]  * 

(65) 

r\/2 

ni-i  ~  \  1 

(66) 

J  -T|/2  1J 

fn/2 

mij  -  3  Tij  * 

(67) 

-n/2 


The  edge  of  the  spherical  shell  at  6  =  eo  is  considered  to 


be  elastically  supported.  Hence  the  boundary  condition  at  8=0. 


can  be  expressed  as 


mv 

j  ail 

a12 

al3 

G14 

ul(90) 

< 

;  a2l 

i 

a22 

a23 

a24 

«  <e0) 

=  1 

I  |  a 

j  a31 

a32 

a33 

a34 

MV  j 

a41 

a42 

a43 

544 

nll(  V 

MV 

mll(V 


(68) 


j  m12(V  ! 

L  -* 


where  a. .  are  the  elastic  constants  of  the  support.  At  the 


apex  0  =  all  displacements  and  stresses  must  satisfy  regularity 
conditions.  The  inelastic  deformation  of  the  shell  will  be  de¬ 
termined  by  the  governing  equations,  Eqs.  (37-67),  the  boundary 
condition,  Eq.  (68),  and  the  regularity  conditions  at  e  =  tt. 
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* 

i 

i 


j 


i 


3.  AXISYMMETRICAL  DEFORMATION 
When  the  magnitude  of  the  external  pressure  is  small,  the 
deformation  of  the  shell  is  axisymmetrical  with  respect  to  the 
polar  axis.  In  this  case,  =  cp  -  e12  =  k12  =  r‘12  =  q2  =  t12  =  0 
and  all  physical  quantities  are  independent  of  3.  We  shall  regard 
this  axisymmetrical  deformation  as  fundamental.  In  the  following, 
we  shall  refer  the  quantities  with  a  super  bar  to  the  fundamental 
deformation  of  the  shell.  The  governing  equations  of  the  funda¬ 
mental  deformation  are 


0*1  =  KUi  -  W  x)  , 

(69) 

‘u  =  5  +  nl>:L  +  > 

(70) 

s22  -•  w  +  u^  cct  9  , 

(71) 

K11  =  %i  ' 

(72) 

ic22  =  5>j_  cot  e  , 

(73) 

Su,l  +  "ll  COt  9  "  n22  COt  9  +  h 

-  *1  hi  P^1  =  °* 

(74) 

) 

\ 


\,i  +  \ coC  9  -  (,'n  +  -  !Vu!,i 

-  cot  0  -  p  =  0  > 

*11,1  +  *11  cot  9  -  ;-22 


(75) 

(76) 


For 

shearing 


axisymmetrical  deformation,  the  shearing  stress  t-j2  and 
strain  e^2  vanish.  Hence  Eq.  (52)  can  be  written  as 


...i. . <  ■* 
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e  =  c  t  +  d 

tv  At  A*  tv 

where 


e  = 

"ell 

j  i = 

?n 

3  £  = 

C11 

C12 

r  1 

dij 

At 

e22 

t22 

C21 

C22 

d2| 

L 

.  J 

where  c. .  are  calculated  by  Eqs.  (53),  (54)  and  (56)  by 
=  t...  The  inverse  relation  of  Eq.  (77)  is 


where 


t  =  f  e  -  g  , 

tv  tv  tv  At 

f  =  c”1  and  g  -  f  d. 

^  tv  At  tv 


r  - 

—  — 

r"u 

£11 

n22 

:  €  = 

At 

€22 

in 

*11 

i 

CM 

CM 

IS 

.  ! 

*22 

By  Eqs,  (51),  (66),  (67),  (79)  and  (80),  we  have 

•  ♦ 

S  =  I  e  •  J  , 


where  I  is  a  4  x  4  matrix  and  J  is  a  4  x  1  matrix. 

tv  At 

elements  in  I  and  J  are 


(77) 

(78) 

setting 

(79) 

(80) 

(81) 

The 
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^q/2 

p  q/2 

f12  d^  I22 

f-Xi/2 

I11  ~ 

flld£  *  I12  “ 

I21  “ 

.\ 

=  \  f22d^  ^ 

-n/2 

-T|/2 

-r|/2 

r  "/2 

rn/2 

H 

M 

CO 

II 

*31  = 

'  \  *11  CdC» 

X14  = 

I23  = 

I32  =  ^l  = 

\  ^]2  2 

-n/2 

-r)/2 

rh/2 

wn/2 

0 

ph/2  2j 

*24  = 

I42  * 

■  \  *22cdc  - 

I33  = 

\  fn^  dC>  V" 

X43"J  f  12^ 

-<l/2 

-q/2 

-n/2 

r 

9  < 

rV2 

rn/2 

ph/ 2 

M 

•P» 

II 

f22^  J1  =  , 

\  9xdC>  d2 

=  \  g2d^ 

J3  =  \  2giCdC.. 

-n/2 

-q/2 

-n/2 

-n/2 

J4  - 

(V2 

g2  cdc  * 

(82) 

-n/2 

where 

and 

g.  are  the  elements 

in  f 

rv 

and  £  . 

The  inverse 

relation  of 

Eq. 

(81)  is 

%  __ 

e  =  K  S  +  L  ,  (83) 

**  ^  ***  ** 

where 

K  =  I"1  ,  L=KJ,  (84) 

By  using  Eq.  (83),  the  incremental  form  of  the  field  equations 
Eqs.  (69-76)  can  be  written  as  the  following  matrix  equation?: 

U  =  AU  +  BV  +  C  (85) 

#v  «v  **  **  ** 

and 

DV-EU+F“0,  (86) 

0V  *»*  tV  **  tv 

where  the  prime  represents  the  differentiation  with  respect  to  9, 

T  T 

S'S5!5  "ll  %  “ill  l  "  t"22  £22J,  (S7) 
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and  the  matrices  A,  B,  C,  D,  E  and  F  are  given  in  Appendix  A. 
After  eliminating  V  from  Eq.  (85)  and  (86),  we  obtain 


vfoere 


and 


U*=  M  U  +  R, 

(88) 

M  =  A  +  B  D"1 

E  , 

(89) 

**  ^  r%f 

R  =  C  -  B  D"1 

*0  **  t*  ** 

Ep 

(90) 

The  boundary  conditions  at  0  =  0Q  can  be  expressed  as 


®X(0o> 

an  a12  al3  j  ! 

nll(so> 

Veo5 

a21  a21  a23  j 

I 

q^Bo)  .  (91) 

i 

5  (Sq5 

i 

a31  a32  a33  |  ! 

—  —i 

!  mu<e0> 

1 -  - * 

At  the  apex,  the  boundary  conditions  are 

^(ir)  =  S(ir)  =  q(7T)  =  0.  (92) 


In  the  following,  we  shall  solve  Eqs.  (88),  (91)  and  (92)  by  a 
modified  Euler’s  method  with  a  finite  difference  scheme. 

Let  n  be  an  integer  and  A  =  (tt-0  )/n  ,  0,  =  0  +  (i-l)A, 

HI  v  M  1  U 


(i  =  1,2 . .  n^  +  1).  U.  =  U(0.),  M.  =  M(0.),  ^  =  R(ei).  Eq.  (88) 

can  be  expressed  as 

s  i+i  -  k + 1  +  &<-i  k+ 1 +  & +  w-  •— v 

(93) 

Let  I.  be  a  unity  matrix.  Put 
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V.  =  I0-  ,  M.  , 

<vi  2  —i 


W  A 

h  ■  io+  2  Si  ’ 

2d.  =  "2  {S±  +  Si+P  ' 
Si  -  ill  1  &  • 


2d.  =  iSil  2d.  - 

Equation  (93)  can  be  written  as 


Jk+1  “  £i  Ik  +  Si1 


(i  =  l,2,...,nm)  (99) 


Hence 


Si  =  Sl  h  +  &  • 


2i“ioand  gj.  =  2 


(100) 


(101) 


By  substitution,  we  obtain  the  following  recurrence  formula  for 
determination  of  a.  and  3.  . 

(%iL 


~i+l 


P.  a. 

!*£L  «*1 


(102) 


&L+1 


=  &L  Si  +  %. 


(103) 


At  the  apex,  we  have 


U_ 

m+1 


■  +  \+i 


(104) 


Hence,  by  the  boundary  conditions,  Eqs.  (91)  and  (92),  we  have 
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all  nll(1)  +  a12  ^1(1)  +  a13 

a21  \l(1)  +  a22  ql^  +  a23  ”ll^ 
•  •  • 
a31  ”ll^  +  a32  ^l^  +  a33 

nn(  1) 

%(!) 

ijjU) 


Ho+i 

m 


0 

0 


5  (nm+l) 


Su'V11 


"ii'V15 


(105) 

The  first,  second  and  fifth  equations  in  Eq.  (104)  are  the  simul- 

•  •  • 

taneous  equations  for  determination  of  0^(1),  q^(l)  and  m11(l)  . 

•  • 

After  U,  is  calculated,  U.(i  =  2,3,...,  n  +1)  can  be  determined 

by  Eq.  (100).  To  avoid  the  singularities  involved  ir.  ci.°  matrices 

A,  B,  and  E  at  the  apex,  we  consider  that  in  the  neighborhood  of 

the  apex,  the  shell  is  under  a  uniform  contraction.  Hence,  at  the 
•  •  •  • 

apex,  we  have  o^ir)  =  h22(7r),  m^T)  =  m22(7r)  and 

I  I  It  J.  t  X  !  it 

cp^(i r)  =  w  (7r)  —  n^ir)  =  q^(T)  =  m^(Tr)  =  0.  From  these  conditions. 


we  can  find  the  matrices  M(tt)  and  R(7r). 

rw*  rw 

Our  calculation  is  based  on  an  iterative  procedure,  i.e.  in 
each  incremental  step,  we  assume  a  condition  of  either  (i)  loading 

a  • 

or  (ii)  unloading  or  reloading  and  find  the  increments  Ui  and 
for  each  spatial  station.  Then,  we  check  the  assumed  condition. 
This  iterative  procedure  continues  until  the  calculated  condition 
agrees  with  the  assumed  condition  at  each  station  of  the  shell. 
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It  is  found  that  the  matrix  involved  in  the  evaluation  of  n^(l), 

•  • 

q^Cl)  and  nL^(l)  is  nearly  singular.  To  improve  the  accuracy  of 
our  computation*  the  computing  program  is  written  in  double  pre¬ 
cision.  The  inversion  of  a  nearly  singular  matrix  is  carried  out 
by  an  iterative  procedure  based  on  the  residual  derived  from  the 
approximate  solution  in  each  iterative  step. 

4.  ASYMMETRICAL  BIFURCATION 

When  the  external  pressure  is  large,  the  deformation  of  the 
shell  may  bifurcate  from  the  fundamental  axisymmetrical  configuration 
to  a  shape  with  an  asymmetrical  mode.  In  the  following,  we  shall 
employ  Shanley’s  concept  of  inelastic  bifurcation  under  increasing 
load  [1]  to  analyze  the  asymmetrical  bifurcation  of  the  deep  spher¬ 
ical  shell.  Shanley’s  concept  was  later  elaborated  on  by  Hill  [3,4] 
based  on  the  argument  of  uniqueness  of  solution  incorporated  with 
the  comparison  theorems  of  linear  and  nonlinear  solids.  Hill’s 
conclusion  can  be  stated  as  follows:  During  the  incremental  process 
of  deformation  of  a  solid  body,  if  the  constitutive  relation  is 
definite  and  independent  of  the  situations  of  loading  and  unloading 
(or  reloading),  the  relations  between  stress  rates  and  strain  rates 
are  linear.  In  this  case,  the  material  is  referred  to  as  a  linear 
solid  and  the  potential  function  is  a  quadrate  of  the  strain  rate 
tensor  with  coefficients  dependent  on  the  state  of  stress.  On  the 
other  hand,  if  the  constitutive  relations  are  dependent  on  the  situ¬ 
ation  of  loading  and  unloading,  then  the  form  of  the  potential 
function  would  depend  on  the  strain  rate  tensor  and  hence  the  con- 
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stitutive  relations  are  nonlinear.  In  this  case,  the  material  is 
referred  to  as  a  nonlinear  solid.  According  to  these  definitions 
of  linear  and  nonlinear  solids  it  is  seen  that  the  potential  function 
WN  for  a  nonlinear  solid  would  coincide  with  the  quadratic  potential 
function  for  the  corresponding  linear  solid  in  a  part  of  the 
strain  rate  space.  In  the  region  where  4  W^,  we  have 
for  most  materials.  It  can  be  shown  by  the  comparison  theorem  that 
if  the  difference  of  potentials  W^-  WL  is  a  convex  function  of  the 
strain  rate  tensor  at  every  point  of  th§  solid  body,  then  the  eigen¬ 
state  in  the  bifurcation  of  a  nonlinear  solid  is  reached  at  a 
higher  bifurcation  load  in  comparison  with  the  bifurcation  of  a  cor¬ 
responding  linear  solid.  Hence,  the  constitutive  relations  used 
for  the  analysis  of  bifurcation  with  the  smallest  initial  load  must 
be  linear  in  the  sense  just  described,  i.e.  the  constitutive  relations 
should  be  definite  and  independent  of  the  mode  of  bifurcation.  From 
the  computational  point  of  view,  the  moduli  used  for  the  incremental 
deformation  due  to  bifurcation  should  be  determined  first  and  the 
bifurcation  analysis  is  simply  a  check  as  to  whether  the  bifurcated 
deformation  is  possible.  Such  a  process  is  similar  to  what  we  use 
for  the  bifurcation  analysis  of  elastic  solids.  Note  that,  under 
an  increasing  load,  it  would  be  proper  to  use  the  constitutive  re¬ 
lations  obtained  from  the  incremental  fundamental  deformation  of  the 
solid  body  for  the  analysis  of  the  critical  condition  of  bifurcation. 

In  the  bifurcation  analysis,  let  us  denote  the  additional  dis¬ 
placement  components  introduced  by  bifurcation  by  u^  and  w,  the 
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rotations  by  <p.  and  m  ,  membrane  strains  by  e..,  bending  strains 
by  • ,  membrane  forces  by  ru  .,  transverse  shears  by  q^,  moment 


resultants  by  nu^,  etc.  We  have 

U1  =  “l  +  “i*  u2  =  -2  '  w  =  ^  +  w 
tpj_  =  i>i  +  cg^  =  S>2  ’  9  =  £  > 


ell  ~  ell  +  ^11  *  e22  “  e22  +  —22  '  e12  “  ^12  ' 


^11  —  ^11  ^  ^  ^22  k22  *-~22  9  ^12  —12  9 


(106) 


nll  ~  nll  +  -11*  n22  "  n22  +  -22  '  n12  "  -12  ' 


ql  ~  ql  +  ql  '  ^2  “  ^2 


m 


11  “  mll  +  -l.V  m22  **  ™22  +  -22"  m12  “  -12 


Substituting  Eq.  (106)  into  Eqs.  (3  7-50)  and  neglecting  the  second 
order  terms  of  those  quantities  due  to  bifurcation,  we  obtain  the 
governing  tinear  equations  for  the  eigenstate  of  bifurcation  where 
the  the  critical  load  is  involved  implicitly  in  the  nonlinear  elastic 
pre-bifurcation  deformation.  In  the  following,  we  shall  assume  that 
g  is  small  and  drop  all  terms  involving  eg  in  the  field  equations.  To 
eliminate  the  P-dependence  in  the  field  equations,  let 

u^( 9,0)  =  OjCe^in  n0,  u2(e,0)  =  u2(e)  cos  n0,  w(9,0)  =  w(e)sin  n3  , 

^(0,0)  =  op1( 9 )  sin  n0,  a2(e,0)  =  cp2(e)cos  n0 


* 


^(0,?)  =  en(0)sin  nP,  e22(0^)  =  e22(e)sin  np, 

e^O, P)  =  e12(e)  cos  np  , 

KiiCe,?)  =  1^(0  )sin  np,  jc22  (0,P)  =  k22(0)  sin  nP, 

JCi2(0,P)  =  K12(0)cos  np  ,  (107) 

n1]L(0,P)  =  nu(e)sin  nP,  n22(Q,&)  -  n22(e)sin  np, 

n12(9,P)  =  n12(e)  cos  nP  , 

Qx  (0jP)  =  <^(0)3  in  np,  qg  (0,P)  =  ^(Q)  cos  np, 

mu<0, P)  =  m11(e)sin  np,  m22(0,P)  =  m22(e)sin  nP, 

migCOjP)  =  1^2(9)003  np, 

where  n  is  an  integer.  The  resulting  field  equations  in  our 


eigenvalue  problem  are 

Cf^  =  X(ux  -  w^)  ,  (108) 

cpg  =  X(u2  -  nw  esc  0),  (1C9) 

eH  “  w  +  ulfl  +  j-  0^  >  (110) 

e22  =  w  +  ui  cot  9  "  nu2  csc  9  3  (111) 

e12  =  "2  1  +  nUl  CSC  9  ”  u2  cot  6  +  T  ^2^  (112) 

*11  =  'h,!  ’  (n3) 

k22  -  -ncpg  csc  0  +  cot  0  ,  (114) 

K12  =  1  ^2,1  +  ^  CSC  6  "  cp2COt  0)  '  (115) 

nll  1  +  oot  0  -  nn12  csc  0-n22  cot  9  +  q1 


“  (cpj_  ^11  ^j_)”  ( HO ) 
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nn22  +  2n^2  dot'  9  +  \2  1  sin  9  +  <l2sin  9’  ^22tp2  +  ^n12)r>in  e 

\  r  * 

-  pepsin  8=0,  (117) 

q^l  +  q-jCOt  e  -  nc*2  cse  9-  (n^  +n22)  -(cpjn^i  +  «llilcp1 

+  *1,1  nll  +  "n^.,15  "  (Vll+  "ll*L>cot  6 

“f  (co^n.^2+  n22®2^  ^  esc  8  =  0, 


mll  1  +  mllCOt  ®  ~imi2  csc  ®“m22  cot®“  \  %  =  0  3 


(118) 

(119) 


nm22  csc  9  +  2m12  cot  9  +  m12ji  ~  j 


i  %2  =  0 . 


(120) 

Since  the  prebifurcation  deformation  is  axisymmetrical,  we  can 


set  t.  .  =  r..  for  i,j  4  1,2  and  t19  =  0  in  Eqs.  (55),  (57),  (58) 
•-*3  X3 

and  (61)  vfaen  we  determine  the  material  constants.  Hence,  if  we 


denote 

- 

”  *1 

f 

i 

nll 

ell  j 

i 

n„„ 

22 

e22 

S  = 

and  e  = 

mll 

K11 

f 

22 

K22 

_  — 

~  - 

(121) 


then 


e  =  X  S  , 

*4  .«W 


(122) 


where  K  is  given  by  Eq,  (84).  The  rest  of  the  constitutive 
relations  are 


and 


'12 


12 


=  *u n 


12 


—  3^22 


19  * 


(123) 


(124) 
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where 


*11  = 


=  1  ±  v 

\T) 


*22  ~ 


-  +  v) 

X  T)3 


(125) 


Hence  Eqs.  (123)  and  (124)  are  elastic  relations.  Let 


and 


U  -  [cp-j.  w  u2  q1  n.^  m12] 


V  -  [c?2  ^  n22  m22  n12^  * 


T 


(126) 


(127) 


Our  governing  equations  can  be  written  as  the  following  equations: 


*  T  *  *  * 

U  =  A  U  +  B  V 


(128) 


and 


*  *  *  * 

D  V  +  E  U  =  0. 


(129) 


*  *  *  * 

The  matrices  A  ,  B  <  D  and  E  are  given  in  Appendix  B.  After 
* 

elimination  of  V  from  Eqs.  (128)  and  (129),  we  have 


*  ?  *  * 

U  =  M  U  , 


(1?0) 


whore 


*  *  *  * 
M  =  A  -  B  D  •lE 


(131) 


Equation  (130)  can  be  expressed  as  the  following  difference  equation: 


where 


*  *  * 

U. =  P.  U.  , 

**1+1  «*a  **i 


*  A  *  -1  A  * 

£d=  <£o-*Sl+i> 


(132) 


(133) 


Thus,  if  we  write 
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we  have 


and 


11?  =  a.  u*  , 
~1  *0.  ~1  ' 


* 

2x  =  I 


*0 


•k  *  * 

a. =  P.  a. 

<*L+1  <*<1  i+OL. 


(134) 


(135) 

(136) 


* 

Hence  all  0L(i  =  1,2,  ...rH-1)  can  be  calculated.  Intuitively,  we 
know  that  when  the  configuration  of  the  spherical  shell  is  deep  enough 
and  the  supporting  edge  is  sufficiently  rigid,  bifurcation  would  not 
be  governed  by  the  case  n  =  1.  The  boundary  conditions  for  cases 
n  ;/  1  can  be  written  as 


U  = 


and 


all  nll<60>  +  a12%(60)  +  Vll^o)  +  a14ml2(9o> 

a31  nll(  V  +  a32^1^9o^  'h  a33mll^0o^  +  a34m12^9o^ 
a21  nll^9o^  +  a22Cll^90^  +  a23mll^9o^  +  a24m12^8o^ 
a41  nll<9o>  +  a42ql^9o^  +  a43mll^  90  ^  +  a44ml2^9o^ 

nn(  0o5 

mll(0o) 

ml2(e0) 

+1  =  [0  o  0  0  q^TT)  nn(ir) 

m 


(137) 


(138) 


&  * 

JL,  +1  -  i  4.1  U, 

«4r^+l  k^ti^+1  «*1 


Since 


> 


(139) 
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with  a  substitution  of  Eqs.  (137)  and  (138)  into  Eq.  (139),  a 
characteristic  equation  can  be  obtained  by  setting  the  determinant 
derived  from  the  coefficients  of  nj_](0o)j  Qi(90)>  mi2^9o^ 

in  the  resulting  first  four  equations  zero.  The  critical  pressure 
for  bifurcation  is  then  determined  by  the  resulting  characteristic 
equation.  In  our  computation,  at  each  incremental  step,  we  check 
the  value  of  the  characteristic  determinant.  The  critical  pressure 
is  determined  from  the  determinant  versus  pressure  curve  for  each 
value  of  n. 

5.  RESULTS  AND  DISCUSSIONS 

In  our  numerical  computations,  we  deal  with  a  deep  spherical 
shell  whose  geometry  is  similar  to  the  hull  of  the  ALVIN  vehicle 
used  for  deep  ocean  exploration.  The  radius  of  the  middle  surface 
of  the  shell  is  40  inches  and  the  thickness  is  hQ  =  1.93  inches.  The 
shell  contains  a  hatch  with  polar  angle  &c  15°.  The  actual  hull 
of  the  ALVIN  vehicle  also  has  three  viewports.  However  in  our 
shell  model,  the  openings  of  viewports  are  neglected  for  simplicity 
of  analysis.  The  shell  is  thickened  in  the  vicinity  of  the  hatch 
opening.  We  assume  a  linear  variation  in  shell  thickness  in  the 
region  15°  s  §  i  24°.  At  the  edge  of  the  hatch  opening,  the  shell 
thickness  is  3.6  inches.  Although  the  computing  program  is  written 
to  include  the  general  case  of  elastic  support  at  the  edge  of  the 
shell,  in  our  computation,  we  only  consider  a  special  case  of  built- 
in  edge  for  which  a^  =  0  (i  =  1,4  and  j  =  1,4)  in  Eq.  (68). 


*5*  ■  1  » 1  "IP' 


W 
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The  material  of  the  shell  is  SK l  2Cb  ITa  -  0.8  Mo  titanium 
alloy  whose  inelastic  mechanical  behavior  can  be  represented  with 
good  accuracy  by  a  Ramberg-Osgood  stress-strain  relation  with 
E  =  16.5  x  lO^psi,  v  -  0.3,  cQ  -  17.  Ox  104  psi,  and  nr  =  50,8.  In  our 
study,  the  effect  of  creep  under  high  external  pressure  is  also  con¬ 
sidered.  Because  of  the  lack  of  the  creep  data  for  titanium  alloy 
under  room  temperature  condition.,  we  set  oc  =  65,4  x  10  psi  and 
m  =  9.9  in  the  power  law  of  steady  creep  fn  our  calculation.  A 
computing  program  is  written  in  FORTRAN  language  for  the  IBM  370 
computer.  The  list  of  the  program  is  shown  in  Appendix  C.  The  time 
interval  used  in  cur  computation  is  so  chosen  that  the  difference  in 
deflection  of  the  shell  is  insignificant  when  we  halve  the  time  inter 
val. 

The  shell  is  submerged  into  the  sea  water  with  a  rate  of  285 

feet/minute  (123.6  psi/min. )  and  then  remains  thereafter  at  a  depth 

of  13,200  feet  where  the  corresponding  external  pressure  acting  on 

the  middle  surface  of  the  shell  is  p  =  0.00742.  When  the  magnitude 

of  this  external  pressure  is  sufficiently  small,  the  deformation  of 

the  shell  is  axi symmetrical.  The  average  normal  displacement,  w  ^ 

is  plotted  against  the  time  in  Fig.  2  where  the  logarithmic  seals 

is  used  for  time.  It  is  found  that  w„  „  increases  as  a  linear 

ave 

function  of  time  for  p  <  0.00742  and  then  remains  at  a  constant 
value  when  p  =  0,00742.  For  this  pressure  history,  the  effect  of 
creep  is  not  noticeable  for  t  <  105  minutes. 
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In  Figs.  3  and  4,  the  tangential  component  of  displacement  in 
the  meridional  direction  uu  and  the  normal  outward  component  of 
displacement  w  are  plottel  as  functions  of  0  for  p  =  0.005194  and 
p  =  0.007420.  It  is  found  that  in  the  region  distant  from  the  edge, 
the  deformation  of  the  shell  is  essentially  the  combination  of  a 
uniform  contraction  and  a  rigid  body  translation  in  the  direction 
of  the  polar  axis.  Hence,  near  the  apex,  u^  and  w  can  be  expressed 
approximately  as 


ux  =  y  sin  0  ,  w  =  5  -  y  cos  8  > 
where  5  and  y  are  constants. 


(140) 


The  membrane  stresses  n^  and  n22  and  the  transverse  shear 
are  plotted  against  0  in  Figs.  6,  7  and  8.  Stress  concentration 
is  found  in  n, .  and  q,  at  the  clamped  edge  of  the  shell.  Due  to 
the  junction  of  the  uniform  shell  thickness  and  the  variable  shell 
thickness  in  the  region  close  to  the  edge,  the  variation  in 
is  not  smooth  in  the  vicinity  of  @  -  24°.  When  the  value  of  0  is 
sufficiently  large,  both  and  n22  approach  a  value  corresponding 
to  the  membrane  stress  in  a  uniformly  contracted  shell  and  n22 
-  p/2  and  the  value  of  approaches  zero.  Tire  variation  in 
and  ir>22  are  snown  in  Figs.  B  and  9.  Again,  it  is  seen  that  both 
r^il  and  approach  zero  when  8  increases.  From  Figs,  5-9,  it 
is  found  that  in  the  region  65°  s  0  s  180°,  the  shell  is  under  a 


uniform  contraction.  The  effect  of  the  donned  edge  ?Ls  only 
restricted  to  the  region  8  <  65°. 
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In  order  to  study  the  effect  of  creep,  the  spherical  hull  is 
submerged  to  the  depth  16,000  feet  and  18,000feet  below  the  sea  level. 

The  corresponding  pressures  are  p  =  0,00900  and  0,01010,  For  these 
pressures,  the  v  versus  t  curves  are  shown  in  Fig,  2,  It  is 
found  that  the  increment  in  the  average  normal  deflection  due  to  the 
steady  creep  is  small.  When  p=  0,01010,  it  is  nearly  3%  for  t  -  10 5 
minutes , 

Next,  we  consider  the  case  without  creep  in  which  the  pressure  can 
increase  indefinitely.  The  pressure  versus  average  normal  displace- 

! 

merit  curve  for  axisymmetrical  deformation  is  shown  in  Fig.  10.  At 
p  =  p  =  0.020701,  equivalent  to  a  depth  of  submergence  of  35,000  feet, 

V 

the  p  versus  wave  curve  reaches  a  maximum  point  which  defines  the 
ultimate  load  for  axisymmetrical  collapse  of  the  shell.  Note  that 
this  ultimate  load  for  axisymmetrical  collapse  of  the  shell  in  the 

| 

inelastic  range  is  much  smaller  than  the  classical  elastic  buckling 
load  of  a  complete  spherical  shell  which  is  (P)cxassicai  ~  2/[3(l-v  )] 

=  1.225,  The  critical  load  for  yielding  of  a  complete  spherical  shell 
is  (PJyiejd-jjjg  =  2td  -  0.0206  which  is  very  close  to  our  pc.  In 

I 

order  to  determine  whether  the  shell  may  bifurcate  before  the  ul¬ 
timate  load  is  reached,  we  employ  the  method  described  in  Section  4 

i 

for  the  analysis  of  asymmetrical  bifurcation  with  increasing  load, 

\ 

It  is  found  that  asymmetrical  bifurcation  can  occur  at  p  =  = 

\ 

0.01916  and  n  =  2.  The  critical  depth  of  submergence  for  bifurca¬ 
tion  is  approximately  33,000  feet.  As  a  result  of  asymmetrical  bi¬ 
furcation,  the  actual  load-carrying  capacity  will  be  smaller  than 

1 

the  one  found  based  on  the  axisymmetrical  deformation.  In  order 


I 


3* _ _ _ _  •  »■>-••  »  - 
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to  determine  the  actual  load- carrying- capacity  of  the  shell,  pest 
bifurcation  analysis  must  be  considered.  Note  that  the  difference 
between  pc  and  p^  is  small  and  we  could  anticipate  that  tne 
actual  load-carrying-capacity  in  the  post- bifurcation  stage  would 
lie  somevAiere  between  these  two  values. 

6.  CONCLUDING  REMARKS 

The  following  conclusions  can  be  drawn  about  the  buckling  of 
deep  spherical  shells  in  the  inelastic  range* 

(1)  When  the  magnitude  of  external  pressure  is  sufficiently 
small;  the  deformation  of  the  shell  is  axisymmetrical.  Asymmetrical 
bifurcation  occurs  when  the  external  pressure  reaches  a  certain 
critical  value. 

(2)  In  the  prebifurcation  stage,  the  deformation  of  the  shell 
in  the  region  distant  from  the  edge  of  the  shell  is  a  uniform  con¬ 
traction  superposed  on  a  rigid  body  translation  in  the  direction  of 
the  polar  axis.  Stress  concentration  occurs  in  the  vicinity  of  rhe 
edge. 

(3)  Under  room  temperature  condition  ,  the  effect  of  creep  is 
usually  insignificant  within  the  limit  of  time  of  operation. 

(4)  For  the  deep  spherical  shell  fabricated  of  titanium  alloy, 
asymmetrical  bifurcation  occurs  when  the  inelastic  deformation  of 
the  shell  becomes  prominent.  The  actual  load-carrying-capacity  of 
the  shell  after  bifurcation  is  smaller  than  the  limit  load  determined 

by  the  axisymmetrical  deformation  of  the  shell.  However,  the  difference 
is  found  to  be  small. 
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APPENDIX  A:  MATRICES  A,  B,  C,  D,  E,  and  F 

w  M  'y  <v 


The  matrices  A,  B,  C,  D.  E  and  F  are  given  as  follows : 
^  ^  ^ 


0 

0 

0 

*31 

0 

CO 

>r 

-cp/X 

0 

1 

•x 

K11 

0 

X13 

-i/x 

1 

0 

0 

0 

0 

Rn  +  p 

0 

0 

5^-cot  9 

-1  * 

0 

T 

"ll+  "ll(^l+  cot6) 

0 

0 

1+^l+K31Flll 

-ii^-COtQ 

K33nll 

0 

0 

0 

0 

1/X 

-cot  9 

_K32 

*34 

• 

i 

K12 

K14 

1 

i 

’ 

0 

0 

j 

I 

i 

cot  0 

0 

t 

\ 

> 

l+J1cot9+K32n11 

! 

V 

lll\ 

i 

L° 

cot 

9  J 

f1®! 


V|-l  V-  IMX-wW  -  > 
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^  A  ^  * 

APPENDIX  B:  MATRICES  A  ,  B  ,  D  and  E 


*  *  * 

The  nonvanishing  elements  in  A  ,  B  ,  D  and 


=  -  1/X  , 


*n  - 


=  -n  esc  9  , 


=  cot  9  , 


r*H_l  +  ^("ll  +  p)  +  nll  cot  6  ' 
“CDi  “  COt  0  > 


-  1  +  ®1,1  +  ^1  +  ' 


K33  nll  ' 


nll  +  P  * 


=  ^  -  cot  0  , 
-  1/X  , 


?w 
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d*  = 
a13 

2n[x(R11  +  K^)  - 

•  Kg2]  esc  0 

d*  = 
a14 

-  2*  <*34 

+  ^22  ' 

■  X^)  esc  9 

d*  = 

15 

-2XK11(©1 

-  cot  9)-  2X^11^1 

d*  = 
21 

1  , 

d*  = 
33 

*22  1 

d*  = 
34 

*24  , 

d*  = 
a4l 

n  esc  0  , 

d*  = 
a43 

*42  ' 

d*  = 

44 

*44  ' 

d|l  =  crv^  -i-  cot  0  , 

=  -2X  Ru  , 

s*i  =  n  esc  0(2  cot  Q  -  3^)  , 

2 

e12  =  n  CSC  8  <2  COt  6  “  X)  > 

e*3  =  ~2n  x  esc  6  cot  9  , 

el4  =  *  ' 

e|6  =  2n  esc  8  (X  -  K^)  , 

e*?  =  2rt  esc  9  (X  l^g  -  Kg3>  , 

e18  =  2  *22  (®!  -  cot  9)  +  2*22,1 

e*  =  n  x  esc  e  , 
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APPENDIX  C:  LIST  OF  COMPUTER  PROGRAM 

THIS  PROGRAM  IS  IN  DOUBLE  PRECISION. THIS  IS  INOICATED  BY 
IMPLICIT  R£AL*8IA-H,0-Z) 

THE  FOLLOWING  IS  THE  INTEGRATION  SUBROUTINE  USED  IN  THIS  PROGRAM 
SUBROUTINE  INT  (V.DEL.N.VAL ) 

NI  *  N-l 
VAL  *  0.0 

DO  201  I  »  l»Nl  ,2 

VAL  »  VAL  ♦  V(I)  ♦  4.0  *  VII  ♦  1)  ♦  VIH-2) 

201  CONTINUE 

VAL*VAL  *  DEL  /  3.0 

RETURN  . . . . .  - 

END  _ _ _ 


EN0 

C  THE  FOLLOWING  IS  THE  DIFFERENTIATION  SUBROUTINE 
DIMENSION  V(51 ) »VSI51 ) 


CD  *  0.5  /OEL 
DO  208  I  *  2.N 
VS ( 1 1  *  OD*?V( 
208  CONTINUE 


i+u  -  vu-in 


vsm  »  od  *  i  -3.o  *  vm  ♦  4.0  ♦  vizi  -  vi3n 

VS  IN  ♦  1)  »  DD  *13.0  ♦  VIN  ♦  i)  -  4.0  *  VIN)  ♦  VIN  -111 

RETURN 

END 

C  THE  MATRIX  ADDITION  SUBROUTINE  FOLLOWS 
ciiPtimiTifc  nMinrw a.a.b.m.m) 


SUBROUTINE  OMAQDIA.B.R.N.M) 
IMPLICIT  REAL *8 1 A-H.p— 2 ) 
DIMENSION  All) .Bill »Rt  1) 


NM  *  N*M 
DO  10  I  >  1,NM 

io  Rin  -  aid  ♦  bid 

RETURN 

END 

THE  MATRIX  MULTIPLICATION  SUBROUTINE  FOLLOWS 
SUBROUTINE  DKPRDIA.B.R.N.M.L) 

IMPLICIT  REAL*8IA-H,0-Z) 

DIMENSION  All) ,8(1). Rll) 

IR  3  0 
IK  «  -M 
DO  10  K  *  l.L 
IK  3  IK  ♦  M 
00  10  J  3  1,N 
IR  3  IR  ♦  1 
JI  *  J-N 
IB  3  IK 
RIIR)  -  0 
DO  10  I  «  l «M 
JI  3  JI  ♦  N 
IB  3  IB  ♦  1 

10  RIIR)  3  R(IR)  +  A( JI )*BI IB) 

RETURN 


TI51 ) .PHI 51  1,0(51  I |Wl5i  ),XNTHI51  ), 
1  ) tXMPHI 51  )» TUTHI 51*21) »TUPHI 51*21) 


DIMENSION  02(5if21)fGI5l»2l)tGSl5i»2l)»CI2f2)|LAI2)*MA|2St0DI2)t 
lGAI2).FJl2,2.2l).GJ12f2l>,Xiil2U,Xl|l2l  ,X|2[21),XI3121).X14  21) 
OIMENSION  X2M21)  tX3lliniX34l2l],X44(|lJ.YUh)  fY2l2l).Y3l21). 

1  Y4  ? 21 ) *THPH (51*21) t PHTH! 51* 21 )  #  TUYftDI 31^21 »  .TUPHDI 51 ,21 ) »CC ( 6 ) 
OIMENSION  Al(4,4),AJI4)1LB?4)1M6l42»ALI4)fAM!6,6),BMI6|2)5CMI6) , 
lXNTHP(5i),PHP{51),gM|2.|),EMl2*6).FM<2)  »DE<2.6)?BDg  6,6) 

lUMD(6i5n  ,UNotit!LO(f3?jMo!S!)!vW&l6?.8N|«6),DAL(3)|GADI3)|ALfC(|) 

lllllliilllltlllilillFi, 


1UMDI6 


6l HENSl6N''5cif  if  Ja ) » LAS!  5  >  * 

dWIA!  f  SAsWiViht  ud?  4?i  irssilf  A? 


WWSI8.B) 
4) , MGAS2 i 


noon 


THE  FOLLOWING  PARAMETERS  TO 
THO  «  HALF  ANGLE  OF  HATCH, 
N  *  NUMBER  OF  INTERVALS  ON 


MERIOll 
TH| 


IN  WHICH  THE  THICKNESS  VARIES  NEAR  T 
THICKNESS  AT  THE  HATCH  EDGE, TO  ■  THE 
BECOMES  CONSTANT, PO  -  VAtUf  QF  PRESS 
NONOIMENIONAL  I 


RATIO,  TUC  _ „ 

TUO*NONDIMENSIONALI ZED  STRf 
SM  ANO  SN  ARE  EXPONENTS  ll 
RESPECT! VELY, XL A  *  RATIO  0 
TD=TIME  INCREMENT,TF«  FINA! 
THE  FOLLOWING  ARE  PROGRAM  Ci 
I  NO  •  0  FOR  NLO  .CREEP  ■I’jfiL  _ 
0 ,  XNTH ,  XNPH ,  XMTH  ,  XMPH  AfftTlH 
EACH  STATION  ON  THE  MERIOI 
IOO  READ  (5,*)  THQ,M.N,NB,H 
IF  KIP  «i  BIFURCATION  J 
FIRST  VALUE  OF  THE  BIFURCA 
THE  INCREMENT.  N81FUR  DENG 
PASS  PRIOR  TO  CHECKING  FOR 


READ  IN  ARE  DEFINED  HERE 
NUMBER  OF  INTERVALS  THROUGH  THICKNESS, 
AL,  NB  «  NUMBER  OF  INTERVALS 
HATCH, HO  IS  THE  RELATIVE 

0J  PRES5U^H^lTA^H^SI^OIN¥^XNU*SPOIisON,  S 
ED  STRESS  CONSTANT  IN  THE  CREEP  RELATION, 
CONSTANT  IN  THE  RAM8ERG-0SG00D  EQUATION, 
THE  CREEP  AMO^RAMBERS-OSGCGD  EQUATIONS 
SHELL  THICKNESS  TO 


MID-SURFACE  RADIUS, 


THE  NUMBER 


MA] 

TER  NIF 
OF  TIME 


1a)«8£nW’ 

I  ARE  NOT  PRINTED. 

IHO,  IP 

NuT*  Nli  DENOTES 
THE  LAST,  N1NC  IS 
INCREMENTS  WHICH 


_  __  _______  _  _  URCATION 

READ  (5,*)  Ail,A12,A13,A14«A2l, A22, A23, A24, A31 * A32,A33,A34,A41, 
1  A42 , A43 » A44 
IF  (M)  101, 500,101, 

101  WRITE  (6,8)  THO,MtN,NB,HO,TO,PO.XNi 
8  FORMAT  (IX,7HTHETA0»E13.5,IX,2HM«I 
i  3HH0«Ei3s5,iXl3HT0*Ei3.5,iXt3HP0« 

I,lX17HLAMB0A*El3,5,lXt7HQfCTAT|EJ 


«  1X,2HN*I3,  lAtanniD* 
13.5, IX,3HNU*E13.5) 


IX , 3HNB*I 3 


IX, 


,  ,  lAt  inu»nou»*ci3i3|  _ .... 

WRITE  (6,2551 )TF, IN3,|P,KIP,N 
2551  FORMAT ( IX, 3HTF-El3.5,ix,6HlN0 
I  4HNI 1*1 3 » IX , 4HNIF* I3,2X ,7HNB 
WRITE  (6.1449) 

1449  FORMAT  ( 2X .32HELASTIC  SUPPORT  CONSTANTS  A(I,J)) 
WRITE  (6,1550)  AU,A12,A13,A14 
WRITE  (6,1550)  A21,A2|,A?3, A24 
WRITE  (6.1550)  A31.A32.A33.A34 
WRITE (6, *550 )  A41 . A42, A43, A44 
1550  FORMAT  (1X,4E13.5) 

FIND  »  IND 

XLB  *  1.0/XLA 

SF33«XLA/(1.0+XNU) 

XN  *  N 
XM  «  M 

TAX  *  3.141592654  -  THO 

DEL  ■  TAX/XN 

N1  *  N  ♦  1 

N2  »  N+2 

NB1  *  NB  +  1 

Ml  *  M  +  1 

DTT*DT 

FNB*N8 

HH* (H0”1 .0 ) /FN8 
TIM  *  0.0 
P  *  0.0 
PR  *  PO/TO 
PRS*0.01*PR*PR 
WAVE  =  0.0 
NT  *  0 

XNN*T0/DT*0.01 
NN  *  XNN 
TUCM  *  TUC**SM 
SMI  -  SM  -  1.0 
SNI  «  SN  -  1.0 
CQ«3.0*SN/7 .0 
DEL2  *  0.5*DEL 
DO  38  I  «  1,6 
R(I.Nl)  •  0.0 
CM( I ) *0.0 
8ET:iI,l)*0r0 
DO  37  j  *  1,6 
AMU  ,  J)*0.0 
ALPI ( I ,J,i)«0,0 
UNK ( I , J )  •  0.0 

37  CONTINUE 

38  CONTINUE 


:I|.5,IX,3HMS*E13.5,1X,3HNS*E13.5 

I.*  I 

X|4HKIP*I3 ,2X, 


,  .a, 


40 


371 


33 

32 


1500 


1501 

1502 


1506 


00  371  1*1,6 
UNMU,li*1.0 
ALPIif ,J,l)*l.O 
CONTINUE 
DO  32  J*1 ,2 
OF I ( J»Nl )  *  0.0 
DO  33  1*1,6 
BMU,Ji*0.0 

mtlknU'0 

continue 

DO  1502  1*1,8 
DO  1500  J*l, 8 
AMSi!t,J  1*0.0 
ALPIS ( 1 , J,1 )  *  0.0 
UNNMI I , J )  *  0.0 
CONTINUE 
00  1501  J*1 *  5 
BMS ( 1  * J)*O.G 
EMS<J,I)  *  OvO 
CONTINUE 
CONTINUE 
DO  1506  I  *  1,8 
UNN^iJ.l)  *1.0 
Ai>iS(t,l,l)  *  1.0 
CONTINUE 

DEI ( 1 ,4, Nl )  *  1.0 
DEI (2,6,N1 )  *  1.0 
AMU, 3)  *  -1 .0 


-XLB 
1.0 
-1.0 
XLB 
0.0 
1.0 
0.0 
0.0 
0.0 
0.0 
»  -XLB 
*  1.0 
*  -1.0 
*  -1.0 
■  XLB 

*  XLB 

*  XLA 

*  -XLA 

*  -1.0 
1  »N1 


10 

11 

12 


AMO,  1  i  * 

AMO, 2)  * 

AMO, 5)  * 

AMO, 5)  * 

EM(1,1)  * 

EMM, 3)  * 

EMM, 5)  * 

EM(2,2)  - 
EMU, 3)  * 

EMt2,5 )  * 

AMSU,l)  = 

AMS (2  r 3 )  « 

AMS (3, 2)  « 

AMS (6,5!  i 
AMS (7 ,5 )  • 

BMS (8,2)  ’ 

EMSU.4)  . 

EMSI2.4)  * 

EMSO.2)  » 

00  14  I  » 

II  *  I  -1 
XU  *  II 

THU)  *  THO  ♦  XI14DEL 
IF  (I  -  NB1)  10,10,11 
T(U  *  HO  -HH*X 1 1 
GO  TO  12 
Tin  *  no 

DELL!  I )  »  TUI/XM 

AKBI (1,1,1)  *  l.0/(SF33*T( I ) ) 

AKB I M , 2  *  I )  *  0.0 

AkIi ( 1 *2 i { ]  -  12?0/(AKBI(l»itI)4T(I)*T(I)) 

DO  13  J  *  i,Ml 

illll.J)  *  -  0.5*T ( ! )  ♦  IXJ  -  l.Q)*0ELL( I ) 
TUTHU,J)  *  0.0 
TUPHM.J)  *  0.0 
TUTHDl I , J) 

TUPHDM  ,  Jl 
J)  -  0.0 
»i  i  » J)  »  c.o 

;S(!,j)*o.o 

IM|I ,  J] 

INTI  NU* 

3NT I  NUf 

i>0  7*41  1*1 ,N 
SI I*OSIN(THt 1 ) ) 


*  •  w 
8:8 


0.0 


8i?T?c2s{T8Jl! 


7941 


7890 

403 

15 

16 
7749 


17 

16 

19 

7704 

7701 

7702 
7711 

7703 

7706 

240 

411 

278 

20 

21 

22 


535 

536 

23 
2  4 
25 


CTU.WCC0/SU 
00  7890  I  «  1»N! 

PH ( 1 1  »  0.0 
UUi  *  0.0 

xAiHtnvi°5.o 

XMPH(I)  *  0.0 

am  *  o.o 

XMTH(I)  *  0.0 
XMPH(I)  *  0.0 
CONTINUE 

WRITE  <6.15)  TIM.P.WAVE 

FORMAT  (Ix15j(TlME*il3a5,lX*2HP-*E13.34lXa5HHAV£*E13.5> 
IF  < IP )240»240» 16 
WRITE <6 .77691 

OTkifP7;i'i*um’ 

FORMAT  ( 1X.4HU( I ) ) 

WRITE  (6.18)  (U(j),I  *  1  ,N1 ) 

rXW - 

FORMAT  i 1X«4HW( I ) ) 

* l*"11 

FORMAT (  •  PH( I )  *») 

WRITE(6.18)  (PH( I ) » 1*1 »N1 ) 

WRITEl6,770lj 

1  1 .  N1 ) 


41 


('  0(1)  »  *) , 
6, is!  (oin.H 

6,7702) 

I1  XNTHU)  -  •) 


(TIM  -  TF)  411,411,100 
NT-NN)  20v278*2l 


FORMAT 
WRITE  < 6 , 

WRITER, 

FORMAT (<  . _ 

WRITE (6.  18J  <XHTH(l>,!*l.Nl> 

WRITE<6,7711) 

FORMAT <»  XNPH(I)  *  •) 

WRITE(6. 18)  ( XNPH( i ) *  1*1 ,N1 ) 

WRITE (6,7703) 

FORMAT ( *  XMTH(I)  *») 

WRI TE (6 , 18 )  (XMTH(I),I*1,N1) 

WRITE  16,7706) 

FORMAT ( *  XMPH(I)*  *) 

WRITE (6, 18) (XMPH( I),IS1,N1) 

NT  *  NT  ♦  l 
TIM=TIM+DT 
IF 
IF 

DT  *  5.0*0TT 
DP  «  PR*DTT 
GO  TO  22 
DP  *  0.0 
DT«2.0*DT 
PST  *  P 
P  *  P  +  DP 
PDOT  *  DP/DT 
ITER  *  1 

tALL  SLOP  (XNTH,DEL,N,XNTHP) 

CALL  SLOP  iPH,DEL,N,PMP) 

DO  24  I  *  1 »N1 
DO  23  J  *  l  •  f*X 

XJ2(I,J)  *  TU(H(I,J)*(TUTH(I,J) ~TUPH( I  *  J ) )+TURH( I , J )*TUPH( I , J ) 
XJM«1,J)  *  DMAXl(XJM(I,J),XJ2(I,  J) ) 

TUE  «  DSQRT(XJ2( I » J) ) 

F  a  Yyg*4SMl/TUCM 
£TA*CQ*( TUE/TUO ) **SN1 

*HPH(I,J)  *  2,o*tuthTi,j) 

PHTH ( I  *  J )  *  2.5*TUPH(1,J) 

XLAr  *  0.5*F/XLA 
01  !  1 ,  J  S  *  XLAF*THPH( I,J)*FIND 
D2  ( I ,  J  )  *  XLAF*PHTH( I, J)*FINO 
IF  (NT-1)  535,535,536 
GSSU  .  J)«0.0 
GO  TO  23 

GSSU.J)  *  2.25*ETA/XJ2(  I, J) 

CONTINUE 
CONTINUE 
DO  36  1*1, N1 

DELT  *  DELHI 
CO  2d  J 


-  TUPH(l,J) 

-  TUTH ( ! , J  ) 


•  »  A 

.Hi) 
1  *  Ml 


G(  I  * J) 

C  C 1 ,2) 
C  (2 , 1 ) 


J)  =  GS( I • J) 

1)  *  XLB*I l *04G( I  * J I *THPH{ IfJ)*THPH(I»J)/9.0) 

2)  *  X^B*t-XNU  ♦  G(  I , J )*THPH( I » J )*PHTH( I » J ) /9.0 ) 

oh  i  nv2  i  s )  the*&Su8le°pr€C  11 !  8n*h atS  i  I  *  invefis  ion  *  ^ubrout  I NE 


CALL  OMINV|C»2jD£Ttl.A»MA) 

OD ( 2 )  *  D2(l,J) 

CALL  DHPRD(C,DD,GA,2,2,1) 


1,2 
1,2 

J)  *  C ( J1 » J2) 


*  X22U> 

FJ(2,2,J)*XSH  I  ,J) 
X13 ( J ) *XSI ( I , J  ) 


DO  27  J1 
00  26  J2 

26  CON^Uui 

G J ( Jl, J )  *  GAU1) 

27  CONTINUE 

28  CONTINUE 
DO  29  J  *  1  *  HI 
X 1 1 ( J )  *  FJil.l.J) 

XI 1 1 ( 1 , J )  -  Xll(J) 

X12( J)  *  FJ( 1 ,2. J) 

X12 1  ( I  ,  J  )  •  X12IJ) 

X22 ( J)  *  F J ( 2,2, J ) 

X22I C I  ,  J I  •  X22U 

X 1 3  <  J )  * 

Xl4(J)  * 

X24U)  * 

X33 ( J  J  * 

X34 ( J )  *  Xl4{J)*XSJI|»J) 
X44(J)  *  XZ4( J)*XSI( 2, J) 

Yi ( J )  s  GJ(l.J) 

Yini.JI  *  Yl  ( J ) 

Y2(J)  *GJ(2, J) 

Y2I ( !  *  J)  *  Y2JJ) 

Y3  (  J )  »  GJU,J)*XS!U»J) 

Y4(  J)  *  GJ(2,J)*XSIU,J) 

29  fKim  (Xll  ,0£LT ,H,  AI  ( 1, 1 )  ) 

SfViJ?  Vp°!*T’*’*,cl’in 

CALL  I NT  (X22, CELT, H«A1 ( 2,2) I 
CALL  INT  (X13,DELT,M,AI(I,3)) 
A I  ( 3  *  1  )  *  A I  ( 1 ,3  ) 

CALL  INT  (X14. CELT, M, Aid, 4)1 
A I  (2,3)  *  A I (1,4) 

AI  (3,2)«Aia,4) 

AI(4,1)  *  AI (3,2 ) 

CALL  INT  {X24.0ELT,H,A!(2,4) J 
AI (4»2 )*AI (2,4) 

CALL  INT  (X33,06LT,H,AI(3,3) ) 
CALL  INT  ( X34.0ELT ,M,AI(3,4) ) 
AI (4, 31  *  AI (3,4) 

CALL  INT 


CALL 

CALL 

CALL 

CALL 


INT 

INT 

INT 

INT 


AI ( 3,4) 

(X44»0ELT,HjAl 14,4) ) 
(Yl,0ELT,H,A4tl) 
IY2,pELT,M,AJ(2 j 
(Y3,0ELT,H,A4I3 


( Y4»0ELTjH,  AJI4)  i 
CALL  DMINV  UI,4,0ET,LBtHB) 
CALL  DMPRO  (AI*AJ,AL,4,4,1) 
DO  300  II  *  1,4 
ALI(Iitl)  »  AL(Il) 

DO  301  12  *  1,4 
AII(I1,I2,I)  *  AI  1 1 1 , 12) 

301  CONTINUE 
300  CONTINUE 

IFU-Ni)366,367,36 
36?  DO  565  11*1,6 
DO  564  12*1,6 
YHK (J 1 , I  2 ' *0  eO 

564  CONTINUE 

565  CONTINUE 
YHK(2f3)*-1.0 
YHK (2,4) *  A 1 1 1,1! 
YMK(2,6)*AI(I,3 
R(2,Ni )*AL( 1 ) 

GO  TO  566 

366  AH (1,4)  *  A 
AH { 1 ,6 )  *  A 


iil:IS 


42 


AMI2, l )  *  -PH<I)/XLA 
AM(2»4)  =  A I ( 1 # 1 ) 

AM (2 (6 )  «  AI  (1*3) 

COP  »  DCOSIPHl!)) 
AM(4,1} 

AMI4.4) 

SP  «  OSINIPHI I ) ) 
AMI5, 1 )  »  XNTHPI 
1  PST*($P  -  PHU) 


XNTHU1  ♦  PST*COP 

ph( i )  -  crm 


4) 

i\ 

6! 


AM  (5 
AM  ( 5 
AM(  5 
AM  <6 
BMU 
BM( 
BM( 
BM  (2 
BM(4 

me  5 

BM(5 


CM<4 
CMC5 
0M(  1 
DM(  1 
OMi2  _ 

iS!i:ii 

EM(1, 4) 
EMU 

m 

EM(2*6) 
FM( 1 )  - 


1.0  4 


m 


HP| 

CT ( I 
)  *XNTHC 


♦  XNTHU  )*IPH(  1)  ♦  CTCIU  - 
‘  ►  PHI  1 )  *PH<  1 ) 


♦  A I ( 3 • 1 )*XMTHI I ) 


A 
All 


8P 


l.o  ♦  p«m*cTm 

4)*XNT><U) 


♦  AI ( 3t 2 )*XNTH( I ) 


♦  AU3»  •  XNTHt It  ♦  PHI I)*PDQT*$P 


*  -At  i _ 

*  -AH2.3 

*  CTjn  . 

*  -AI  4,1) 
«  -AI 14,3) 
ALI2) 

AL 14  J 


CALL  unruu  ion,ur  ,Duri\,o,i. 
CALL, DMADO I AM^  BDE « YMK ,6,6) 

CMU1)  -  BOFKl  Ill 


rn«wr*c»2< . 

DF , BDFK.6, T , 1 ) 


35 


)0  35  II  * 

%\inU*  c 

>0  203  II  * 


•  infill 
2  «  1,6 
2,1)  »  OE ( 1 1 • 12  ) 


11 

12 


DO  35 
R 
C 

00 

OF! Ill, 

00  202 
DE! Ill, 

202  CONTINUE 

203  CONTINUE 
566  00  40 

DO  39 
OELY-DE 
VMUl.I 
Will, I 

39  CONTINUE 

40  CONTINUE 

CALL  OMINV  IVM,6,0ET,LC,MC) 

00  42  II  »  1,6 

00  41  12  «  1,6 

VMM 1 11,12,1)  *  VMI II, 12) 

41  CONTINUE 

42  CONTINUE 
36  CONTINUE 

DO  45  I 


1,6 
1,6 
L2*YNK( 1 1 
2)  *  UNM( 

*  I  >  *  UNM 


iit I21  -  0 

Hill, 12)  4 


OELY 
DELY 


00  44  II 

44  ^ONtInUE 

45  CONTINUE 
00  50  I 
IA-i+l 
DO  46 
BET(  ID 
zzcii) 

DO  47  12 
ALP  111,12) 


‘fi 


L2*(RI  Kill)  4  R  { 1 1 , 1  4  U) 


1  »N 


n  >  1.6 

l-BETIill.I) 
*  ZMIII,!) 


1.6  . 

•ALPI(U,!2,t) 


43 


44 


) 


WWU1.I2)  *  WM (11,12, I) 

VM (11*12)  «  VMN( 1 1 , 1 2, 1  A ) 

47  CONTINUE 
46  CONTINUE 

Eut 

CALL  OMPRO  ( VW , AL  P »  ANe *6*6*6) 

CALL  ONPRO  ( VW*BET ,VWB,6,6, 1 ) 

DO  54  11*1,6 

BETim.iAj-vwBumvzun 
DO  53  12*1,6 

ALP  I (11*12* I  A) -ANE (11*12) 

53  CONTINUE 

54  CONTINUE 
50  CONTINUE 

X0(1)  *  0.0 
XG(2)  *  0.0 
Xu(  3 )  *  0.5 
IT-1 

GAM( 1«1)-ALPI(1,4,N1) ♦All*ALPI(  1 « 

^AAtl!  2)-ALPI ( 1 ,5»Nl >4A12*ALPI( 1, 

1.3. N1) 

GAN(1,3)«ALPI(1,6,NI)4A1**ALPI(1, 

5aA(2  ,  l)-ALPI(2*4*Nl)  4All*ALPI (  2  * 

1.3.  Nl) 

GAM(2.2)-ALPI(r  5,N1)*A12*ALPI(2, 

1.3,  Nl) 

GAM(2.3)=ALPI(  ,6,N1)*A13*ALPI(2, 

1.3,  Nl) 

GAH(3,1)-ALPI (5,4,N1)+A11*ALPI(5, 
I ,3»N1 ) 

ALPI (5,6,N1)+A13*ALPI(5, 


=  ALPH5,5,N1)+A12*ALPI(5, 


GAM(3, 3) 
l,3,Nl) 

UAL ( 1 ) -BET I ( 1 *N1 ) 

DAL (2 )*BET! (2*N1 ) 
DAL(3)-8ET1(5,N1) 

DO  91  11-1,3 
DO  90  12-1,3 
GBM1 11,12)  *  GAM( 11,12) 

90  CONTINUE 

91  CONTINUE 

CALL  DMlNV  tGAH93.OET1LO.KD) 

CALL  OHPRO  (GAM, DAL, GAO, 3,3, 1 ) 

735  X0(1)  «  X0(1)  -  GAOl 1 ) 

X0I2)  »  X0(2)  -  GAD(2) 

X0( 3 )  *  X0( 3 )  -  GAD( 3 ) 

CALL  DHPPD(GBK,X0»Y0»3,3,1 ) 

DSL  ( 1 )  *  YOU)  ♦  OAL(I) 

DBL (2)  *  Y0(2)  ♦  PAL (2) 

DBL(3)  «  Y0( 3)  ♦  OALCl) 

AAA  *  DMAX1(0ABS(DBL(12) ,DABS(DBL 
IF (AAA-  0,00001)  798,790,799 

799  IT  *  IT  +  I 

I F  < I T  -  10)  732,732,733 

733  WRITE  (6,734)  IT 

734  F0RMAT(1X*3H IT-13) 

732  CALL°OM?RO(GAM,OBL,GA0,3,3,i) 

GO  TO  735 

798  UMD (1,1)  -  A1 i*XO( 1 )  ♦  A12*XQ(2) 
UMD (2,1)  -  A21*X0tl)  ♦  A22*XO  2) 
UMD(3,1)  -  A3l*X0( 1 )  +  A3**XQ( 2 ) 
UMD (4,1) *X0( 1 ) 

UMO ( 5 , 1 } *X0(2 ) 

UMD (6*1 ) *X0( 3 ) 

DO  788  II  •  1,6 
CC(ID-  UHD(  11,1) 

788  CONTINUE 

CO  59  1-2, Nl 
DO  57  11*1,6 
DO  56  12-1,6 

ALP ( 1 1 , 1 2  J-ALPI (11,12,1) 

56  CONTINUE 

57  CONTINUE 


l#Nll  )4A2l*ALPI  ( 1,2,  Nl 
1 iki  >+mz+M.p  Itl ,  2  ,  N 1 
l*fUl*tt3*ALPI(i,2,Nl 
l,Ni )+A21*ALPI ( 2,2,Nl 
1,N1)*A22*ALPI(2,2,N1 
1,N1)4A234ALPI(2,2,N1 
l,Nl )+A21*ALPI ( 5,2, Nl 
l,Nl)4A22*ALP! ( 5,2,N1 
1,N1)+A23*ALPI(5,2,N1 


)+A31*ALPI(l 
)4A32*ALPI(1 
)+A33*ALPI (1 
)+A31*ALPI ( 2 
)+A32*ALPI (2 
)+A33*ALPl (2 
)+A3l*ALPI (5 
)+A32*ALPI (5 
)+A33*ALPl  (5 


(21), DABS ( DBL  <  3 ) ) ) 


♦  A13*X0( 3 ) 

♦  A23*X0(3) 
+  A33*XO(3) 
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CALL  DMPRD  ( ALP,CC, ALPC,6,6, 1 ) 
00  58  f 1*1 *6 

UH0( 1 1 »  I ) *ALPC( 1 1 ) ♦BET  I ( 1 1. 1 ) 

58  CONTINUE 

59  CONTINUE 
DIF*0.0 
OEO  *  0.0 

DO  77  I  *  l »N1 
CK ( I }  *  0.0 
CWK ( I )  *  5.0 
DO  67  12  *  1.6 
UNO ( 121  *  UM0U2«I} 

B2i ii*  t2»xi 

66  CONTINUE 


67  CONTINUE 

11  * 


DO  68 
DF(I1) 


n 


68  CONTINUE 

CALL  DMPRD  <0E,UND.VD,2,6, 1 ) 

$D(  2 )  =  vom  -  DPI  1 ) 

SD<4)  *  V0I2I  -  Or (2) 
void, I)  *  SD(2 ) 
von?, n  *  so <41 
soil)  «  UND(4) 

SD<3)  *  UN0<6) 

DO  70  II  *  1,4 
DO  69  12  *  1.4 

Aim, i2)  *  A:im,i2tn 

69  CONTINUE 

70  CONTINUE 

CALL  DMPRD  ( A I , SQ,A$D,4,4, I ) 

DO  71  II  *  1.4 
ESD(Il)  =  ASD(Ii) 

71  CONTINUE 
DO  ?6  J  *  l, Ml 

py  :  m\i\ :  n\\M\ :  mu 

tuthd< i , j)  *  xniir 

TUPHDU.J)  *_Xf 


♦  AL I < 1 1 «  I ) 


XJ2DU.J)  •  THPH1 J  ,  J)*TyTHO(  I ,  J 1  ♦  J  >*TUPHD<  I , 

IF  (DABS ( XJ2D( I »  JH  -  PRSf  72,72,73 

72  G ( I  * J 1  *  0.0 
GO  TO  75 

73  IF  (  XJ2D!  I  •  J 1 1  72.72,744  ..  _ 

744  IF(XJ2U,J)*DT+XJ2<I,J)-XJM<  I,J1)  Y2*74,74 

74  G(  I ,  J)  *  GSSUyJ) 

75  CKD  -  DA8S(G(I,J  -  G$ C I  * J > ) 


CKK  *  DABStGl I . J; 1 
CWK ( I )  *  OMAXl (CWK( 1 1 ,CKK ) 

C K  <11  *  OMAXl (CK( I ) , CKD) 

GS(I.J)  *  GU.J) 

76  CONTINUE 

DIF  *  DMAXKCM  I). DIF) 

OED  •  DMAX1 (CWK( I J ,DEul 

77  CONTINUE 

IF  (DIFi  81,81,777 
777  IF  (OIF/DED  -  0.001)  81,81,78 

78  ITER  *  ITER  ♦  1 

IF  < ITER  -  20)  25,25,79 

79  WRITE  (6.80)  .  . 

80  FORMAT  { IX , 18HITERAT IONS  DIVERGE ) 
GO  TO  500 

81  I F { KIP )  1081. 1081, j 
1102  IF(NT-N8IFURi  1081, 

1100  NCOUNT*NII-NIN* 

7192  NC0UNT*NC0UNT*NINC 

IF {NCOUNT-NIF >1787,1787, 1081 
1787  SNN*NCOUNT 

00  1101  I»1 »N1 


1102 

,1081,1100 


110 


)-AKBjU,i*n 
}*AKBI<2,2,I> 


AKBC 
AKC  ( 

CONTI 

CALL  SLOP  (AM 
CALL  SLOP  (AK< 

00  7175  I»i,N] 

I F ( I-Nl )  1788,1789,1789 
1789  00  1790  II  *  1,8 


EL,N,AKCP) 


J ) 


DO  1791  12 
YMKS( 1 1  *  T  ? 1 
1791  CONTINUE 
1790  CONTINUE 
GO  TO  7194 
1788  AMS { l »6) 
AHS(1,7) 

AMS (3  , 1 ) 
AMS'  3,6) 

AMS (3,7) 

AMS (4*3) 

AMS (4,4) 

AMS (5*1) 

AMS (5*5) 

AMS (5,6) 

AMS (5*71 
AMS(6, 1) 

AMS ( 6  »6 ) 
AMS(7  ,7) 

AMS (7,8) 


«  1.8 

•  0.0 
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Ain3.i»t) 
AII(3.3,I) 
“PHI  I ) *XL8 
Aiiiiti.n 
Aim.S.n 
-SNN*csm 
cTm 


XNTHP(.l|«;Pt(i I )*(XNTH( I ) ♦PST  > ♦XNTHf I l*CT ( 1 ) 
I.0+PHPm4pH(n*PH«  I )+AI 1(3, 1, I  )*XNTH(  I ) 

feuirTHm 


ph  m -crm 
-CTm 
SNN*CS(n 


Mi’.h:f;rCTin 

00  1503  J JJ*1 ,5 
DMS(1 1 1, JJJ)«6.0 
1503  CONTINUE 
1505  CONTINUE 

DMSU.ll*-l-0+  PHP(n*2.0*PH(I)*CT(I)  ♦  (  AK8I(  2.  i ,  !  )“XLA*AKBI  1 1 , 1 . 1 
1 ) ) * C XNPH (I) +PST ) *2.0 

DMS(1,2)  *  2.0*U^B*  AKBS ( 2,  2.  !  ) “AKB 11 1.2.1)  “I  AKBI  <  2, 1, 1 )“ 
lDMS?!!3lli“2?0*SNN*CSm*tAKBM2tl»n-IAKBm.l,n*Aini,2,n  )*XLA 

lDMsllU?,*>-2.0*SNN*CS(  I  »*iAXBM2,2*n-XLA*AKBm.2.I)+Am3.4,n” 

1  All  (1,4,1 )*XLA) 

OMS  (1*5)  *  2.0*(  (CT(I)+PH(I) )*AKBI(2»  1*1  )“AKBP(  I  )*XLA*  , 

1 ( AKB  £(1*1*1) *XLA-AKBI (2*1.1) ) *( -PH(  I  )  *2 .0*CT  C 1 1 ) -AKB I(l.l.I)*CT(I) 

2  *XLA ) 

*  1.0 

*  All(2,2,n 

=  All (2*4.1) 

*  SNN*CS ( I j 
«  AI 1(4. 2.1) 

=  AI i (4,4, I ) 

*  CTUJ+PHd) 

*  2;0*(AKBI  (2,1.  D-AKBH  1,1,1  )*XLA) 

»  AI  (3,2, 1 ) 

*  AI  (3,4,1) 

*  AII(l, 2, if 
»  AIKl.  4, 1) 

*  -XLB*PHl I ) 


OMS (2,1 
OMS (3, 3 
DMS(3,4 
OMS (4,1 
DMS(4,3 
OMS (4,4 
DMS (5,1 
0MS(5,5 
BMS (1,3 
BMS ( 1 ,4 
BHS {5 1 3 
BMS(3,4 
BMS (4,1 
BM5(4,5 
BMS (5,1 
BMS (5,2 
BMS (5,3 
BKS(5,4 
BMS(6,3 
8MS(6,5 
BMS (7*4 
BMS(8,4 
6MS (1,1 
EMS (1,2 
EHS ( 1 » 3 
EMS (1,6 
EMS (1,7 
EMS (1*8 

1  -mctIi 

EMS (2*2 
EMS( 3, 3 
EMS (3 ,4 
EMS (3, 6 
EMS(3,7 
EMS (4,1 
EMS (4,6 
EMS ( 4  *7 
EMS  (  5 , 1 
EMS ( 5  *2 
EMS  ( 5 , 3 
EMS (5  *4 


2  »0*AKBI (1,1,1 ) 

-$NN*XNPH( I l*C$( I ) 

SNN*CS(II  _ 

1.0+PH1 I )*CT III *XNTH( I) PA 11(3, 2,1) 
AII(3,4,I)*XNTH( I) 

CT(n 

8MSC5.2) 

BMS (6,3 ) 

“BMS ( 5,2 ) 

SNN*CSi 1 )*( 

SNN*CS ( I ) ♦ 


«-2.0*SNN* 


*CT(  n-3.0*PH(I)> 


)*CT( 


=  2 


•  0*SNN*CS( I) *( AI I  ( 1 , 
.0*SNN*CS(I)»(Aini, 


n*CT(I)“l»0)4XLA 

)*XLA 

1,1 )*XLA“AI I ( 3, 1 , 1 ) ) 
3,II*XLA-AIH3,3,n  ) 


«  2»0*( AKCP( I )-( AKBI ( 2, 2, I )-AKBI ( i,2, I )*XLA )*CT( I )*2.0 
^PH(I)  )*AKBI(2,2,n-CTTI)*AKBI(l,2,l)*XLA) 


SNN*CS ( I ) *XLA 
-CT( 1) 

ara?i!i. 

All (2,3,1 ) 

sillk:  1 1 

“2.0*SNN*CS( I) 
-SNN*CS( I)*CT(!I*XLA 
-EMS (5,1) *XLA 
EMST4,1)*XLA 


7176 

7177 
7194 


EHSCS.S)  *  2.0*{AKBI{2.2*n-AKBm,2»n*XLA) 
CALL  DMINVCDMS*5,DET,LAS,MAS) 

CALL  DMPRD{DHS»EMS»DE5»5»5,8 ) 

CALL  DNPRD(BMS,DES,BDE$,8,5,8) 

00  7177  IJ1  *  1,8 
DO  7176  IJ2  «  1,8 

YHKSC I Jl , IJ2 )  «  AMSC I Jl, ! J2)-80ES( 1J1, I J2) 

CONTINUE 

CONTINUE  . 

00  7172  I Jl  *  1,8 
00  7171  IJ2  ■  1,8 
DELYS  *  DEL2*YMKS( I Jl *  I J2> 

VMS( I Jl , I J2)  »  UNNKC I Jl,  Li2)  DELIS 

NMS{IJi,IJ2,I)  «  UNNM(IJl,?J2>  *  DELYS 
CONTINUE 


NHSC I Jl, I J2, I )  « 

7171  CONTINUE 

7172  CONTINUE 

CALL  DMINV(VMS#8| 
00  7174  IJ1  «  1,8 
00  7173  IJ2  »  1,8 
VMMSMJ1,IJ2,I)  * 

7173  CONTINUE 

7174  C0N7INUE 

7175  CONTINUE 

00  7180  I  «  1 «N 
IA  »  141 

DO  9177  IJ1  *  1,8 
00  9176  I J2  *  1,8 


DET *  LCS»MC$) 


VNS(IJ1,1J2) 


ALPS ( I Jl *  I J2 )  *  ALPISI ljlt I J2, 1 ) 
WHS ( I Jl , I J2 )  *  WMS{ I Jl, ! JZ» I ) 

VMS ( I Jl , I J2 1  *  VMMSI 1 Jl , I J2, IA) 

9176  CONTINUE 

9177  CONTINUE 

CALL  0MPR0(VHS,HHS,VMS,8,8,8) 

CALL  DHPRDC VW$ , ALPS , ANES • 8,8,8) 

DO  7179  I Jl  *1.8 
00  7178  I J2  *  1.8 
ALPIS(1J1,IJ2,IA}  *  ANESdJl*  IJ2  ) 

7178  CONTINUE 

7179  CONTINUE 

7180  CONTINUE 


GAMS(l.l)  *  ALPIS(1,6»N1J+ALPISC 1,1*N1 )*A11+ALPISI l,2,Ni 
l  ALPlSU,3,Nl)*A2WALPfSU,4.Nl)*A41  .  „  . 

GAMS (1,2)  *  ALPIS(i,5,Nll+ALPISU,l,Nl)*A12+ALPIS(l,2,Nl 
1  ALPlS(l,3,Nl»*A22*ALPIS(l,4tNl}*A42 


GAMSC 1.3)  -  ALPISIJL,7,NlJtALPiS4iaxNl 
ALP  I S  Cl  ,  3 ,  Nil  ALP  IS  11 . 4  .  N1 )  **$I 

MS ( 1 ,4 )  *  ALPIS ( 1 ,8«Nl J+ALPISC 1,1, N1 
ALPl  SCI, 3, N1)*A24+ ALPISI  1,4,N1}*M4 
MSC2.1)  *  ALP1S(2,6,N1]+ALPIS(2*I*N1 
ALPIS<2,3.N1)*A21*ALPI$C2,4,N1)*A41 


GAMS (l 


)*A13*ALPISU,2,N1 
)*A14*ALPIS( 1,2,N1 


GAMSC2, 1 )  *  ALP1SI2,6,N1]+ALPIS(2*I*N1>*A11+ALPISC2,2»N1 
1  ALPI§(2,3,N1J  *A21*ALPI S ( 2  »4*N1 )*A4i 

GAMS  (  2,2 )  *  ALPI S  C2 , 5>,N1  l+ALPISC  2*  1  ,N1  )*A12>ALP I  SC  2, 2,N1 
l  ALPI$(2»3»N1) ♦A22+ALP I  SC  2,4, Ni )*A42 

GAMS  12,3)  *  ALPISiSL,  7 jNL) tALPlSC  2 ,1  ,Fil  )*A13+ALP ISC  2, 2,N1 


7183 


7747 

7746 


l  ALPlS(2»3»Nl) ♦A22+ALP I S( 2,4, Ni )»A42 

GAMS  12,3)  *  ALPISi 2, 7jNL) tALPlSC  2,1  «FJ1  )*A13+ALP ISC  2, 2,N1 
1  ALP  I SC2,3»Ni)*AZ34ALPl5C2,4,NI )*«43 
GAMSC2.4)  «  ALPIS (2«8,NiT+ALPiSC  2, l, HI )*A14+AvPI SC  2, 2* N1 
1  ALPlS(2,3,Nl)4A24+ALP|S(2,4tNl)*A44 

GAMSC 3, l )  -  ALPISI3,6,N1 IfALPlSC  3, 1, N1 )*A1 i+ALPI SC  3, 2»N1 
1  ALPI$(3,3,Nl)*A2l+ALPI5i5,4,Nl>*A4l  .  „  . 

GAMS  C  3  *  2  >  *  ALPISC3,5*Nl]4ALPIS(3,ltNl)*Al2+ALPISC3,2,Nl 
l  ALPISC3*3,NI) ♦A22+ALPI § ( 3»4,N1]*A42 

iG,s^?§n.;.ik^iik{!m^UiHbSu*Al34AlPISI5’2*N 

GAMS ( 3*4 )  *  ALPIS^a.NlH-ALPlSM.l.Nl^AU+ALPlSO^Nl 
1  ALPISC3,3,NlJ#A24+AtPISC3,4.Nl)*A44 
GAMSC4.U  *  ALPI ST4,6jW1 )+ALPISC 4, ItNl )*A1 1+ALPI S(4, 2,N1 
1  ALPlS(4f3,Nn*A21  +  ALPlS(4,4|Nl)*A4l 
GAMS (4 ,2 )  *  ALPIS(4,5,N1 ) +ALP ISC  4, 1, N1 )*A12+ALPI $ ( 4, 2 ,N1 
1  ALPIS{4,3,Nl)*A22*ALPf$C4,4.Nn*A42  ,  „  , 

&WAU*AlPIS'Al  ,n 

GAMS (4,4 )  »  ALP I SC 4, 8 »Nl )*ALPlS( 4*1, N1 )*A14+ALPI SC 4, 2,N1 
1  ALPlSC4,3,Nl)*A24+AlPlSi4,4.Nl}*A44 
THE  MATRIX  GAMS  IS  MULTIPLIED  BY  10**“10 


)*A31  + 
)*A32<- 
)*A33+ 
)*A34* 
)*A31+ 
)*A32+ 
)*A33< 
)*A34+ 
)*A31* 
)*A32+ 
)*A33^ 
)*A34+ 
)*A3l* 
)*A32+ 
)*A33* 
)*A34+ 


00  7746 
DO  7747 
GAMSC 11, 
CONTlNUf 
CONTINUE 


10**- 10 


11-1,4 
! 2*1 ,4 

12)  -  G AMS  ill, 12) *0,0000000001 
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uni  n* io*njittuc I 

(6,7195)  NCOUNT ,DET 1 , PST 
TI*  NCOUNT**  ,  I3»  *  DE 


;  CALCULATION  OF  THE 
CALL  DMINVIGAHS 
WRITE 

7195  FORMAT 

GO  TO  7192 
1081  00  83  I  *1 »N1 

UU)  *  U(  1  >  ♦  UNO 
win  *  win  ♦  ukdu» 

PH(I)  *  PHI  l  >  ♦  UM0(Un«0T 
XNTH(I)  *  XNTHIIJ  ♦  UHD(4,I>*DT 
xnph(I)  *  xNPHin  +  voni,n*oT 
QU)  -  Oil)  ♦  UM0I5  — 

XMTHII 
XMPHI I 

Ml 


DETERMINATE  OF  GAMS* DENOTED  8Y 
4*DET1 *HGAS1*NGAS2 T 


T1  *  * , E13.5, • 


DET1 


PST 


m 


82 

83 


DO  82  J  * 
TUTH(I,J) 
TUPHI I  * J) 
CONTINUE 
CONTINUE 


*  ' 

♦  UMDI5.n*0^' 

iMitiHis? 


TUTHC I  * J) 
TUPHI I, J) 


♦  TUTHDI I* J>*0T 

♦  TUPHDl 1  * J 1*DT 


C  CALCULATION  OF  AVERAGE  RADIAL 
CALL  INT<W.DEL,N,WINTT 
WAVE  *  WINT/TAX 
GO  TO  403 
500  STOP 
END 

/* 

//GO. SYSIN  DD  * 


DEFORMATION*  WAVE 


*  .E15.7 ) 


UN L- 73 -7 


- 


X 


3 


Fig.  1  The  geometry  of  the  shell. 
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